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Abstract 



We study the transport properties of a NSN structure with an insulating 
barrier at each NS interface. Coherent quasiparticle scattering is assumed 
and self-consistency is implemented exactly to guarantee local charge conser- 
vation. The presence of a finite condensate flow has a greater influence on 
the transport properties than either the gap depression near the interfaces 
or the coherent nature of scattering. We find that a nonzero phase gradi- 
ent causes a shift towards lower voltages of the first peak in the differential 
conductance and a global enhancement of current. At low currents, we ob- 
tain gap profiles near the interfaces that are consistent with the criteria for 
boundary conditions employed in macroscopic descriptions. The existence of 
coherent multiple scattering gives rise to a rich structure of resonances that 

is smoothed out for long superconductors. 

PACS numbers: 74.40.-hk, 74.50.-Fr, 74.80.Fp, 74.90.-Fn 
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I. INTRODUCTION 



Motivated by the recent development of mesoscopic super conduct ivityjll the study of 
transport in normal-superconductor structures has been the object of renewed attention. 
Microscopic studies are conventionally performed within the framework of the Bogoliubov - 
de Gennes (BdG) equations:!! 
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where [wnj^n] and are the wave function and energy of quasiparticle n. Hq is the one- 
electron Hamiltonian referred to the Fermi energy and A is the pair potential. The BdG 
equations are based on a mean field Hamiltonian which does not commute with the particle 
number operator .1 It has been notecllH§ that current conservation is only guaranteed if the 
BdG equations are solved within a self-consistent scheme, i.e., one in which the gap function 
A(r) is required to satisfy the condition 

A = gY,UnV:il-2fn), (2) 

n 

g being the electron-phonon coupling constant and /„ the occupation probability. 

The requirement of current conservation can be satisfied only if the condensate carries 
a finite amount of currentJBi and for this a nonzero phase gradient is needed: ^ 0, 
with A = lAle*"^. At distances from the scattering region much bigger than the coherence 
length, the superfluid velocity Vg = TiS/(p/2m and the gap amplitude acquire uniform val- 
ues. Based on this result, a model was developed in Ref.i to compute the current- voltage 
characteristics for several normal-superconductor structures with the simplifying assump- 
tions of asymptotic self-consistency (current is globally conserved) and incoherent multiple 
scattering (if more than one interface is present). A rich transport behavior was unveiled 
of which the most salient features are a shift towards lower voltages of the first peak in the 
differential conductance and an enhancement of current at not very low voltages caused by 
the greater availability of current carrying scattering channels. Different voltage regimes 
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appear in the presence of a finite condensate flow, since this creates a distortion in the 
quasiparticle dispersion relation within the superconductor which makes the various quasi- 
particle channels open at different voltages. A regime was found in which quasiparticles 
enter and leave the superconductor only through Andreev transmission and in which cur- 
rent is thus insensitive to the presence of impurities within the superconductor.! It was also 
shown that a state of anisotropic gapless superconductivity can develop for high enough 
voltages if the temperature is low .01 Effects due to a finite superfluid velocity are important 
when the supercurrent is comparable to the bulk critical current. This may occur already at 
voltages of order Ao/e if the interfaces are transmissive and the superconductor is at most 
as wide as the semi-infinite normal leads, i.e., if there is no geometrical current dilution in 
the superconductor.li0 The self-consistent gap near a NS interface with zero current was 
already studied by McMillanS and, more recently, by BruderEl in the context of anisotropic 
superconductivity. 

The purpose of this paper is to go beyond the approximations introduced in the model 
of Ref.i and study the nonlinear transport of a NSN structure within a picture in which 
quasiparticles undergo multiple coherent scattering and current is locally conserved (self- 
consistency is implemented exactly). In particular, we wish to know which of the transport 
properties found in Ref.i survive in a more accurate description. Like in previous scattering 
studies of transport through NS structures,lilii we assume that an insulating barrier may exist 
at the interface contributing HS{x) to the one-electron potential Hq. The dimensionless 
parameter Z = mH/%^kF is a measure of the barrier scattering strength, since the one- 
electron reflection probability is Z'^/{1 + Z"^). We are interested in NSN structures which 
contain two insulating barriers -one at each interface- of strengths Zi and Z2. By performing 
a systematic study of the dependence of transport properties on the values of Zi and Z2, as 
well as on the length of the superconductor, we extend the work of Martin and Lambertlll 
where a similar set of assumptions was employed. 

In an asymmetric NSN structure, the distribution of the potential drop between the two 
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interfaces has to be determined self-consistently, like in the asymptotic model of Refi. Here 
we determine in addition the exact profile of the pair potential A(a;). Consistently with 
previous work^ we find that the position of the first peak in the differential conductance 
(FPDC) is shifted to lower voltages, because of the existence of a nonzero phase gradient. 
This effect is more important than the local depression of A(x) or the coherent nature of 
quasiparticle scattering. Due to the coherent multiple scattering, a structure of resonances 
develops which gives rise to an oscillatory behavior of dl/dV as a function of voltage and 
length (see Figs. 4 and 5). The results obtained from the model of incoherent scattering 
and asymptotic self-consistency are recovered qualitatively for structures where the super- 
conductor is much longer than its coherence length. 



II. THE MODEL 

We wish to solve the BdG equations (1) self-consistently within the context of a one- 
dimensional model. The resulting gap (2) will in general be space dependent. At a given 
energy, a scattering problem has to be solved which is defined by Eq. (1) with Hq describing 
a NSN structure with insulating barriers at the interfaces. Quasiparticle scattering states 
are labeled by the incoming channel a. Thus, index a indicates the lead (left or right) and 
the type (electron or hole) of the initial quasiparticle state. We can exploit the fact that only 
wave vectors in a small interval around ±kp are important, and linearize the one-electron 
Hamiltonian around the two Fermi points. In momentum representation, 

-— = —{±kF + kY Ef ± — = -Ef± ^k. 3 

2m 2m m m 

Therefore, instead of Hq = {Ti^ /2m){—(P' / dx^) — Ef^ we can write 

fi^kp d ^ , 

Ho-±t -—-2Ef, 4 

m dx 

in the vicinity of =F^f- For the numerical resolution, the BdG equations are discretized within 
the superconductor S, so that only a finite set of positions Xi = iL/N (with i = 0,1, ...,N) 
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are considered, with a uniform spacing S — L/N. Within the interval [xi — S/2, Xi + d/2\ the 
solution of Eq. (1) can be written as a sum of quasiparticle plane waves, Z]/3 ^/3Qe*'^^*'^[M^, v^] 
which are solutions for the "uniform" gap A{xi). This is a safe procedure if 5 is chosen much 
smaller than the superconducting coherence length = ^vf/ttAq, which is the typical length 
scale for gap variations. At each point and energy, there are four possible quasiparticle 
states labeled by /3. The components at i + 1 are obtained in terms of the components 
at i through a 4 x 4 transfer matrix that is determined by wave function matching at 
X — {xi + Xi+i)/2. When a linearization scheme is adopted, the ^kp branches do not couple 
within the superconductor (only Andreev reflection or normal transmission can occur) and 
the transfer matrix at each point is box-diagonalized into two 2x2 matrices. Right at the 
interfaces {x = and x — L) normal reflection and thus branch mixing occurs, and in those 
points the transfer matrix is also determined by standard matching techniques. The global 
transfer matrix is obtained by compounding the individual transfer matrices. 

Once the scattering quasiparticles have been calculated, the gap is determined from Eq. 
(2). Due to normal reflection at the interfaces, a quasiparticle n has in general components 
from both ^kp branches. This gives rise to irrelevant oscillations in A(x) on the scale of the 
Fermi wavelength, which we smooth out by introducing only intra-branch interference terms 
in the gap equation (2). Once self-consistency has been achieved, the resulting description 
does conserve charge and the electric current must be uniform. 

Self-consistency is implemented as follows. The input parameter is the difference in 
chemical potentials between the quasiparticles coming from the left and the right normal 
leads. The desired output is the resulting current. The pair potential profile A.{x) and the 
chemical potential at the superconductor ns are treated on the same footing, as parameters 
that are progressively adjusted to achieve convergence. In the first scattering calculation, 
some guess values for the A{x) profile and /is are introduced. From the global transfer 
matrix we obtain the quasiparticle scattering states, and from them the current in the 
normal leads. For the following iteration, the gap at each point is obtained from Eq. (2) 
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and the superconductor chemical potential is adjusted to make the currents in both leads 
equal. Self-consistency is achieved when A(x) and iis converge and the current in the two 
leads is the same within tolerable errors. 

The main simplifications that we have introduced are: (i) Use of a one- dimensional model, 
which we expect to yield qualitatively correct physics in structures with planar interfaces 
where mode mixing is not important, ii) The voltage drop takes place entirely at the 
interfaces and the chemical potential within the superconductor is uniform, (iii) There is 
no proximity effect {g{x) = in the normal metal). A similar set of approximations were 
employed in Ref.,0 with the difference that no linearization procedure was implemented 
there. Relaxation of (ii) and (iii) may yield interesting physics, but here we want to focus 
on the effect of a nonuniform, complex A(a;). 



III. RESULTS 

There are two main ways in which local current conservation can affect the transport 
properties of a NSN structure. One is the induction of a nonzero phase gradient, which has 
already been mentioned. The second effect has to do with the local depression of the gap 
near the interface, which is a direct consequence of self-consistency that already occurs at 
Vs = 0. From the work of Refs.0&0 we have some understanding of the effect of fiow and 
now we wish the explore the effect of gap depression at the interface^ In Fig. 1 we plot for 
comparison the Andreev reflection probability as a function of the quasiparticle energy for 
gap profiles with a step-like structure (solid lines) and with a local, smooth depression shown 
in the inset (dotted lines). The profiles have been obtained by truncating (and matching 
with a uniform gap) the self-consistent gap of a NSN structure with zero current. Thus, 
there is no fiow effect in any of the two cases. A similar study has been performed by 
van Son et al.El introducing phenomenological gap profiles with a proximity effect. The 
influence on scattering of the local gap reduction is small in the cases of low and high Z, 
but it is appreciable for intermediate Z. At energies E ^ Aq, Andreev reflection (AR) is 



6 



enhanced for Z = 0.5 because it takes place in three distinct (albeit coherent) steps: the 
electron is first normally transmitted, then it is Andreev reflected as it finds a higher gap, 
and finally the emerging hole is normally transmitted. AR at the smooth gap increase takes 
place with probability very close to unity. On the other hand, the two normal transmission 
(NT) processes are more probable than direct AR (since this requires the simultaneous 
transmission of two particles), which is the only possibility when E < A{x = 0). Thus, by 
being decomposed into three processes of higher probability, AR takes place more easily. 
The effect of gap depression is very small at Z = 0.1 because direct AR already has a high 
probability. For high Z, the gap depression is small and thus it cannot have important 
consequences. 

In Figs. 2 and 3, we plot the gap amplitude |A| as a function of position and voltage for 
two asymmetric NSN structures: (Zi, Z2) = (0.1, 0.5) (Fig. 2) and (2, 0.5) (Fig. 3). Here, Z^ 
and Z2 are the barrier parameters at x = and x = L, respectively. In both cases, L = 8^o- 
At low voltages, | A| is strengthened if the barrier is thick, while it is more strongly depressed 
near a transmissive interface. This is in qualitative agreement with the Ginzburg-Landau 
result for the boundary conditions describing the interface between a superconductor and 
an insulator or a normal metahli At a transmissive NS interface, the order parameter must 
be zero, while at the interface with an insulator (here represented by a high Z interface), it 
is the derivative of the order parameter what must vanish. We see againS that the physics 
obtained from a self-consistent description at zero temperature is qualitatively similar to that 
which is derived from the Ginzburg-Landau approximation near the critical temperature. 
As the voltage increases, it drops mostly at the more reflecting interface, where a greater 
gap reduction takes place because of the stronger presence of quasiparticles. In fact. Fig. 
2 shows the rather counterintuitive feature that |A| can even become smaller at the more 
reflecting interface, in marked contrast with the zero voltage behavior. 

Upper Fig. 4 shows the differential conductance dl/dV as a function of voltage for 
several structures as obtained from different types of calculations. Thick and solid lines 
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have been obtained with exact and asymptotic (see Ref.S) self-consistency, respectively. The 
other three lines have been obtained from uniform phase calculations, in the spirit of Ref.lll: 
The dotted (dashed) line corresponds to incoherent (coherent) scattering, while the long- 
dashed line results from a coherent scattering calculation performed with the self-consistent 
gap profile of the zero current case: A{x,V = 0). In all cases, Z2 = 0.5, L = 8C,o, and 
T = 2 K have been taken, with superconductor material parameters corresponding to those 
of Pb (Tc = 7.2 K). The lower part of Fig. 4 shows the relevant energies of the problem at 
the midpoint x = L/2 = 4^o corresponding to the exactly self-consistent calculation: |A| 
(dotted), A± = |A| ± hqvp (dashed), and the voltage difference at each interface (solid). 

The first peak in the differential conductance (FPDC) occurs when eV reaches A_, 
and its position is very sensitive to the existence of finite superconducting flow. This is 
particularly clear in the case of the symmetric {Zi = 0.5) NSN structure. The lowering 
of the FPDC is less marked in the Zi = 0.1 case because the asymmetric voltage drop 
tends to mask the effect.! Finally, for Zi = 2, the FPDC is barely shifted, because of a 
global reduction of current caused by poor transmitivitty. As compared with calculations 
based on asymptotic self-consistency, implementation of local self-consistency and coherent 
quasiparticle scattering displaces the FPDC position only slightly, and tends to increase the 
current. 

The existence of coherent multiple reflection at the interfaces gives rise to a structure 
of resonances that is clearly appreciated in the plot of dl/dV vs. voltage. The oscilla- 
tory behavior of the differential conductance predicted by coherent scattering calculations 
contrast with the less structured curves obtained from the assumption of incoherent scat- 
tering. These resonances are very efficient at carrying current because they are made of 
quasi-bound states in which the quasiparticles are mostly Andreev reflected at the inter- 
faces. Thus, for instance, a right-moving electron is Andreev reflected as a left-moving hole, 
which in turn is reflected as a right-moving electron, and repetition of the process results in 
a strong net electron current to the right. These oscillations are very similar in nature to 
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the Tomasch oscillations occurring in tunnel junctions,0 and also occur in other transport 
contexts involving coherent Andreev reflection.lll 

The period of the oscillations in dl/dV is strongly sensitive to the total superconductor 
length L. Actually, for a superconductor with constant A = Aq, one would expect to have 
resonances at energies 
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which is obtained from requiring — = {2n — l)7r/L [k^ are the two possible solutions 
to e{k) = E in the branch, e{k) being the quasiparticle dispersion relation). In our 
case, Eq. (5) is not exactly fulfilled because of the effect of f ^ 7^ on the pair potential. 
Nevertheless, Fig. 5 shows that the separation between peaks decreases with increasing L, 
with a first peak that is not very sensitive to the length (see inset). The spacing between 
peaks depends also on the values of Z (not shown) because of the effect that these have on 
the effective gap amplitude. 



IV. CONCLUSIONS 

We have performed a self-consistent scattering calculation of nonlinear transport through 
asymmetric NSN structures. By introducing local current conservation and coherent scatter- 
ing of quasiparticles, we have gone beyond the model employed in Refs.i'i, where incoherent 
multiple scattering and global current conservation was assumed. Comparison of the results 
obtained from the two different models allows us to identify the most robust physical fea- 
tures which are likely to survive in realistic scenarios. In both models we have obtained a 
lowering of the voltage threshold for quasiparticle transmission (as signaled by the first peak 
in the differential conductance) and a global enhancement of current caused by the increased 
availability of charge-transmitting scattering channels. Both effects are directly related to 
the presence of a nonzero superfluid velocity in the condensate. 

In general, the self-consistent gap does not have a uniform amplitude. Linearization of the 
BdG equations around the two Fermi wave vectors automatically washes out any spurious 
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spatial dependence on the scale of the Fermi wavelength, and only the relevant physics 
occuring at the scale of the coherence length is preserved. The gap amplitude |A(a;)| tends 
to decrease with increasing voltage and is locally depressed in the vicinity of the interface 
with the normal lead. In the presence of a strongly reflecting insulating barrier the depression 
is weak and | A | drops abruptly to zero at the interface (no proximity effect is assumed) . By 
contrast, if the interface is transmissive, the gap amplitude drops smoothly from its bulk 
value to a very small value right at the interface. This result is consistent with the criterion 
for boundary conditions employed in macroscopic, Ginzburg-Landau descriptions.! However, 
we have found that the presence of a nonequilibrium population of quasiparticles may alter 
this picture qualitatively. 

Quasiparticle multiple scattering is responsible for the existence of resonances within the 
superconductor. These quasibound states are very efficient current carriers, since they con- 
sist mostly of electrons that are Andreev reflected at the interfaces, or viceversa. Whenever 
the voltage reaches a threshold level that permits transmission through a new resonance, 
there is a quick rise in the current, hence the nonmonotonous dependence of the differential 
conductance on the applied voltage (see Figs. 4 and 5) . We have found that agreement with 
incoherent scattering calculations tends to improve as the superconductor length becomes 
large compared with the coherence length, because in that case the structure of resonances 
is smoothed out. 

We end this article by pointing out some remaining challenges that should be addressed 
in future scattering studies of current conserving transport. Self-consistency should be ex- 
tended to include the Coulomb interaction term, since voltage variations due to charge pileup 
at the barriers may be comparable to the gap if the interface is appreciably reflecting. Most 
important, it would be desirable to understand how the transport properties discussed here 
and in previous references are affected by the presence of many modes which can mix at 
the interface and by collisions with impurities in the superconductor. Finally, it will be of 
interest to understand how transport is modified by the existence of a proximity effect in 
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the normal lead. 
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FIGURES 




FIG. 1. Andreev reflection at a NS interface with different insulating barriers for a stepwise 
(solid line) and self-consistent (dotted line) pair potential. Insets show the self-consistent pair 
potentials employed in the calculations. 
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FIG. 2. Amplitude of the self-eonsistent pair potential for a NSN structure with Z\ = 0.1 (at 
x = 0) and Z2 = 0.5 (at x = 8^0) as a function of position and applied voltage. A temperature of 
r = 2 K (for Tc = 7.2 K) has been taken. 
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FIG. 3. Same as Fig. 2 for Zi = 2 and Z2 = 0.5. 
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FIG. 4. NSN structure. Z2 = 0.5 and data of Pb for S have been taken in all cases. Temper- 
ature is r = 2 K while Tc = 7.2 K. The upper curves show the differential conductance calculated 
with exact (thick solid line) and asymptotic self-consistency (thin solid). The other three lines 
result from calculations with a zero phase gradient: Dotted (short-dashed) line corresponds to 
incoherent (coherent) scattering with a stepwise pair potential, while the long-dashed line corre- 
sponds to a coherent scattering calculation performed with a zero-current gap. The lower curves 
show the relevant energies (in units of Aq) of the problem in the exact calculation: Voltage drops 
at each interface (solid), magnitude of the order parameter |A| (dotted), and of the A+ and A_ 
thresholds (dashed) evaluated at x = L/2 = 4^o- 



16 





r 



5.0 10.0 15.0 20.0 25.0 30.0 35.0 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 

eV/Ao 

FIG. 5. dl/dV YS. F for the Zi = 0.5 case of Fig. 4 for different superconductor lengths: 
L = 8^0 (soUd), L = 12^0 (long-dashed), L = 16^o (short-dashed), and L = 24^o (dotted). Inset: 
Positions of the various resonance peaks as a function of the superconductor length. 
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